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Abstract 

Numerical results obtained from the analysis of the 
conjugated heat and mass transfer of an Air-CCte mixture 
inside a ventilated cavity in turbulent flow regime have been 
presented. Two configurations (case Dl and D2) of the cavity 
were analyzed, in order to study the thermal behavior and 
the air quality inside the cavity. A contaminant source (C02) 
of 3000 ppm is considered in nine different positions inside 
the cavity. These positions are identified with the name of 
CHI, CH2, ... CH9. The air inlet gap is located on the lower 
side of the vertical wall of the cavity: in the right wall (case 
Dl) and the left wall (case D2). From the results, for the Dl 
configuration, the lowest value of Taverage is obtanied in cases 
where the source of contaminant is located near the heat 
source (CH3, CH6 and CH9). While for the D2 configuration, 
the best performance from the standpoint of thermal and air 
quality corresponds to the cases with a source of contaminant 
in the positions CHI, CH2 and CH3. 
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Nomenclature 

C average'. Average contaminant concentration (CO2) 
covering the inside of cavity. 

Qnkt: Concentration of chemical specie at the mixture 
inlet, Kg-m- 3 o ppm. 

Coulet'. Concentration of chemical specie at the mixture 
outlet, Kg-m- 3 o ppm. 

Csource, Ch: Contaminant source of the chemical specie, Kg-m- 3 
o ppm. 

W: Width of the cavity , m. 
H: Height of the cavity , m. 

Q: Imposed heat flux on the exterior of opaque wall, W 

m- 2 . 

Qconv-ext'. Convection heat flux towards the exterior of the 



cavity, W mr 1 . 

Cjrad-ext'. Radiation heat transfer towards the exterior, W 

m- 2 . 

Re: Reynolds number, Re = (u inlet )(p)(H i ) / ju . 

Taverage'. Average temperature of the mixture (Air-CCh) 
inside the cavity, Q C 6 K 

T inlet'. Temperature of the mixture (Air-CCh) at the 
mixture inlet, Q C 6 K. 

Toulet'. Temperature of the mixture (Air-CCh) at the 
mixture outlet, Q C 6 K. 

UMet: Horizontal velocity at the mixture inlet, m S' 1 . 

Greeks Symbols 

E: Rate of disipation of k, m 2 /s 3 

S t : Overall effectiveness coefficient for temperature 
distribution 

S c : Overall effectiveness coefficient for contaminant 
distribution 

k : Turbulence kinetic energy, m 2 /s 2 . 
Introduction 

Indoor air pollutants derived from a range of sources 
are emitted by the fabric of buildings, and may also be 
a by-product of the activities that are undertaken 
within them. Sources can be broadly classified as being 
associated with the activities of building occupants and 
other biological sources, the combustion of substances 
for heating or fuel, and emissions from building 
materials. The concentration of a pollutant indoors 
depends on the relationship among the volume of air 
contained in the indoor room, the rate of production or 
release of the pollutant, the rate of removal of the 
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pollutant from the air via ventilation, the rate of air 
exchange with the outside atmosphere, and the 
outdoor pollutant concentration (Jones). Ventilation is 
closely related to the comfort level experienced by a 
person. The comfort is a function of multiple variables 
such as: air temperature, average radiant temperature, 
air humidity, air velocity, type of clothing, level of 
activity, etc. Another purpose of ventilation is to 
provide a supply of fresh air to the inside spaces in 
order to dilute contaminants concentration produced 
by people, equipment and materials. Contaminants 
dilution is influenced by the amount and quality of 
exterior air, as well as the way the air is circulated to 
the interior space (Awbi, Allard). Other useful works 
related with natural convection in open cavities have 
been performed by Prakash, Nateghi and Armfield, 
Singh and Venkateshan and Mohamad. This technical 
note evaluates the performance of a room with a 
source of pollutant (CO2) located within it, in order to 
obtain thermal comfort and air quality inside room, the 
traditional way to model ventilated room is thorough 
cavities with different configurations and conditions. 
The ventilation parameters used in the evaluation are 
the overall ventilation effectiveness for temperature 
distribution and contaminants removal. 

Configuration and Mathematical Model 

The physical configuration consists of a two-dimension 
cavity with the air inlet gap on the lower side of the 
right wall (Case Dl) and the left wall (Case D2). For 
both configurations, the air outlet gap is on the right of 
the upper wall. The cavity is ventilated through the 
lower side of the vertical wall with air at temperature 
Tinlet and a C02 concentration Cinlet. The vertical 
right wall is a solid opaque wall subject to a constant 
heat flux falling normal to the (q) surface. Convective 
and radiative losses to the outside of the cavity are 
considered. Finally, the contaminant source (C02) is 
considered in nine different positions inside the cavity. 
These positions are identified with the name of CHI, 
CH2, ... CH9. The first position (CHI) is in the upper 
left side of the cavity at a distance of 1 m from the left 
vertical wall at a height of 2.25 m. The position CH2 is 
horizontally at a distance of 1 m from the position CHI 
with the same height. Position CH3 is at a distance of 1 
m from the right vertical wall at the same height of the 
polluting sources CHI and CH2. Other contaminant 
source positions are in that order downwards at a 
distance of 0.75 m, in equidistant positions, as seen in 
Figure 1. The nine positions of the source of 



contaminant have been considered on the 
configurations Dl and D2. The height and width of the 
cavity is defined as H and W, respectively. The inlet 
and outlet gaps height is H, the inlet air velocity (uinut) 
is a function of the Reynolds number (Re) and the 
conductive wall thickness is Wm. The governing 
equations for mixed convection inside a ventilated 
cavity are mass, momentum, energy and chemical 
specie (CO2) conservation equations, averaged in time. 
These equations systems are not complete due to the 

presence of the Reynolds stress tensor (pu { u . ) in the 

momentum equation, the turbulent heat vector 

(pu t T ) in the energy equation and the turbulent 



mass vector (pu t C ) in the chemicals specie. In the 

Eddy Viscosity Model (EVM), the Reynolds stresses, 
turbulent heat and mass flux vector are approximated 

u ; du . 



as: 



pu i u j = -ju t 



dx. 



2 * 
+ pKO..> 

3 J 



pC u.t =- 



M t dT 



dx i 

and p u \c 



Where 



cr dx i ' 1 Sc t dx i 

<j f is the turbulent Prandtl number and Sc f is the 



turbulent Schmidt number. The turbulent viscosity ([it) 
is related to turbulent kinetic energy (k) and to the 
turbulent kinetic energy dissipation (a) by the 
Kolmorogov-Prandtl empirical expression as: 



PK 



In order to fix the turbulence 



mathematical problem, the turbulent kinetic energy (k) 
and the turbulent kinetic energy dissipation (e) are 
used (Henkes et al., [8]). The boundary conditions for 
the velocity components are on all the solid surfaces 
and they are zeros, at the air inlet gap, they are 
u = u inlet and v inlet = , at the mixture outlet gap 

— = and = , where n is the vector normal to 

dn dn 

the flow direction. The thermal boundary conditions 
are at the air inlet gap T = T Met , at the mixture outlet 

gap dj_ _ q . On the opaque surface in contact with the 

dn 

fluid inside the cavity, from a thermal 
balance, q cond _ muro = q conv . int , and finally for any other 
region on the boundary the walls are adiabatic. 
Boundary conditions for the chemical specie (CO2) are 
at the air inlet gap C = C inlet , and at the mixture outlet 

gap we have = , on the position (CHI, CH2, .CH9) 

dn 

there is a contaminant source that is C = C Source and 
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for any other region on the boundary walls ,it is taken 

dC 

as = • For the right wall, the governing equation 

dn 

is defined as the heat conduction equation. The 
boundary conditions for this model are the lower 
(y = H t ) and upper (y = H ) boundaries are adiabatic; 
and over the boundary in touch with the environment, 
a constant heat flux is imposed with convective and 
radiative losses to the ambient temperature T ext , that 

is to say, q = q cond _ muro + q rad _ ext + q conv . ext • For the 
boundary interacting with the cavity, it is assumed that 
the conductive heat flux through the opaque wall is the 
same as the convective heat flux to the cavity. 
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FIG. 1 PHYSICAL MODEL FOR THE VENTILATED CAVITY 

Characteristic Parameters 

The overall ventilation effectiveness for temperature 
distribution S t , provides a quantitative index related to 

the way in which heat is distributed covering inside 
the cavity, the higher the value is, the more 
homogeneous the temperature distribution is. Awbi 
defined this parameter as _ _ T outlet —T inlet . So do the 

f T —T 

average inlet 

contaminant removal effectiveness S c which provides 



a quantitative index related to the way in which 
pollutants are distributed covering inside the cavity, 
the higher the value of S c is, the more homogeneous 
the pollutants distribution is. This parameter is defined 



as 



s r = 



C -C 

^outlet _ inlet 

C -C 

average inlet 



Numerical Procedure 

The governing equations describing the flow, heat and 
mass transfer were discretized by the finite volume 
technique (Patankar, [9]). Coupling among the mass 
and momentum conservation equations has been 
carried out using the SIMPLEC algorithm defined by 
(Van Doormal and Raithby, [10]). Convective terms 
were approached by the hybrid scheme and diffusive 
terms by the centered scheme. The algebraic equations 
system is solved through the line by line method (LBL) 
with alternating direction implicit scheme (ADI). If the 
values in the mass balance for every control volume as 
well as the residual values of the different equations 
are significantly low, then the overall convergence will 
be achieved (typically 10 8 ). The accuracy of the 
numerical results is verified through numerous tests 
based on the grid size effect. A mesh independence 
study has been carried out with numerical meshes 
varying from 71x61 to 121x111 with 10 nodes 
increments on each axis. The discretization for the 
opaque wall is always 10 nodes on the horizontal 
direction for all the meshes. A maximum deviation of 
1.0 % for the variables between meshes 111x101 and 
121x111 has been observed. Therefore, a 111x101 grid 
was used for all cases considered herein. The results of 
this validation were presented by Serrano [11], with 
satisfactory results. From the comparison of results, 
our numerical code in consideration is sufficiently 
tested. The Reynolds stresses that turbulent heat and 
mass flux vector are approximated as :- 



P u i u j 



dx. 



du. 



PC P u i T 



<j f dx. 



2 * 



P u iC 



Sc t dx. 



(i) 

(2) 
(3) 



Results and Discussions 

For this study, an actual size room cavity 4.0 x 3.0 m 2 
(H x W) has been considered. The inlet and outlet of 
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the mixture air-CCh are H/=0.3 m height, the same size 
as the air conditioning diffusers. The material for the 
opaque wall is construction common brick with the 
thickness Wm=12 cm . Supplied air temperature has 
been set to be 24°C and dniet =340 ppm and the outside 
conditions for the opaque wall are a convection heat 
transfer coefficient (hext) of 6.8 W/m 2 K, which 
corresponds to a wind velocity of 3 m/s and ambient 
temperature of 35°C. The heat flux on this wall had a 
constant value of 750 W/m 2 (ASHRAE, 2005 [12]). The 
Reynolds number is considered for values 10000 and 
20000. The CCh contaminant source considered is 3000 
ppm (maximun values for a room environment). 

Table 1 shows the summary of results for each of the 
locations of the contaminant source (CHI, CH2, CH9 
...) for the Dl configuration with Re = 10000. The 
average temperature and concentration of CO2 in the 
cavity is shown. Similarly, it shows the temperature 
distribution efficiency (et) and the contaminant 
removal effectiveness (£c) inside the cavity. It is 
observed that average temperature values are in the 
range of 24.3 to 24.5°C, corresponding to the position 
of the source of contaminant in CH2, CH3, CH6 and 
CH9, for these positions, and larger values of et can be 
achieved. In addition, the lowest value of temperature 
(24.3°C) has been obtained for the location CH3, with 
the corresponding value et = 77.8%. In general, it can be 
said that from the thermal point of view, the best 
results are obtained for the positions of the 
contaminant source closest to the outflow opening. The 
average concentration is the lowest when the pollution 
source is located in the top positions: CHI, CH2 and 
CH3 (from 382 to 391 ppm). For these positions, too, 
belongs to the highest values of £c (from 40.5 to 58.3%). 
These results indicate that removal of contaminant 
should be at the top, the closest to the position of the 
source of contaminant. Further, it can be seen that the 
locations CH7, CH8 and CH9, corresponding to the 
position of the source in lower part of the cavity, have 
the lowest values of £c (13.1 < £c< 13.4). 

Analogously, Table 2 shows the summary of results for 
each of the locations of the contaminant source (CHI, 
CH2, CH9 ...) for the Dl configuration with Re = 20000. 
It is observed that the lowest value of temperature 
(24.2°C) has been obtained for the location CH9, with 
the corresponding value et = 52.6%. This et result shows 
that when the Reynolds number is increased, the 
temperature homogeneity decreases within the cavity. 
The average concentration is the lowest when the 
pollution source is located in the top positions: CHI, 
CH2 and CH3 (from 403 to 418 ppm). For these 



positions, too, belongs to the highest values of £c (from 
34.6 to 38.1%). These results indicate that removal of 
contaminant should be at the top, the closest to the 
position of the source of contaminant. It is observed as 
well that, this £c result shows that when the Reynolds 
number is increased, the CCh-Concentration 
homogeneity declines within the cavity. 



TABLE 1 RESULTS FOR Dl CONFIGURATION FOR RE = 10000. 



Re=10000 


Csource 

POSITION 


Taverage 

( Q C) 


Caverage 

(ppm) 


et 


£c 


ChI 


24.7 


382 


38.8 


58.3 


Ch2 


24.5 


390 


47.8 


52.3 


Ch3 


24.3 


391 


77.8 


40.5 


Ch4 


25.2 


435 


21.0 


22.8 


Ch5 


24.9 


460 


26.9 


21.5 


Ch6 


24.5 


440 


49.5 


21.4 


Ch7 


25.0 


463 


24.0 


13.4 


Ch8 


25.0 


493 


25.1 


13.6 


Ch9 


24.4 


470 


57.0 


13.1 


TABLE 2 RESULTS FOR Dl CONFIGURATION FOR RE = 20000. 


£e=20000 


Csource 

POSITION 


Taverage 

( Q C) 


Caverage 

(ppm) 


Et 


£c 


ChI 


24.8 


403 


25.4 


34.6 


Ch2 


24.7 


418 


32.7 


373 


Ch3 


24.5 


405 


39.1 


38.1 


Ch4 


24.9 


477 


22.3 


17.2 


Ch5 


24.6 


469 


32.5 


18.2 


Ch6 


24.2 


421 


71.5 


20.8 


Ch7 


25.0 


463 


21.0 


13.7 


Ch8 


24.5 


483 


38.6 


13.5 


Ch9 


24.2 


448 


52.6 


16.1 



Table 3 shows the summary of results for each of the 
locations of the contaminant source (CHI, CH2, CH9 
...) for the Dl configuration with a Reynolds number of 
10000 and 20000. The same ventilation parameters of 
Table 1 and 2 are presented in this table. For Re=10000, 
it is observed that the values of lower average 
temperature are obtained for three locations: CHI, 
CH2 and CH3 corresponding to the cases where the 
pollution source is located at the top of the cavity (~ 
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24.1 to 24.2° C). Of these three positions to the 
pollution source, CHI position has the best 
distribution of temperature (30.5%). The average 
concentration is the lowest when the pollution source 
is located in the top positions: CHI, CH2 and CH3 
(from 384 to 419 ppm). For these positions, too, belongs 
to the highest values of £c (from 49.1 to 54.9%). The 
better distribution of pollutant corresponds to the 
location CH3. Additionally, it can be seen that the 
locations CH7, CH8 and CH9 which correspond to the 
position of the source in lower part of the cavity, have 
the lowest values of £c (12.8 < e c < 13.9). For a value of 
Re = 20000, similar to the results of Re = 10000, the best 
results for ventilation correspond to positions CHI, 
CH2 and CH3. This comparison shows that increasing 
the value of Reynolds number allows an increase in the 
values of et and e c , from which it is concluded for this 
configuration that the homogeneity of temperature and 
concentration of CO2 have been improved. 

TABLE 3 RESULTS FOR D2 CONFIGURATION FOR RE = 10000 AND RE = 
20000. 



&KL0000 



Csource 


Taverage 


Coverage 


£t 


£c 


POSITION 


( Q C) 


(ppm) 






ChI 


24.1 


384 


30.5 


49.1 


Ch2 


24.2 


395 


24.2 


50.6 


Ch3 


24.2 


419 


22.3 


54.9 


Ch4 


24.8 


446 


5.6 


22.6 


Ch5 


24.7 


456 


6.1 


21.5 


Ch6 


24.5 


444 


8.3 


20.9 


Ch7 


25.2 


504 


3.4 


13.9 


Ch8 


25.2 


515 


3.6 


13.6 


Ch9 


24.7 


491 


6.3 


12.8 


#e=20000 


ChI 


24.1 


396 


59.3 


73.1 


Ch2 


24.1 


404 


62.2 


73.0 


Ch3 


24.1 


391 


50.0 


73.7 


Ch4 


24.3 


405 


10.5 


58.2 


Ch5 


24.4 


441 


9.8 


61.2 


Ch6 


24.3 


465 


11.7 


51.4 


Ch7 


25.0 


429 


4.0 


50.9 


Ch8 


25.0 


438 


3.6 


44.8 


Ch9 


24.5 


435 


7.4 


37.4 



www.seipub.org/rame 
Conclusion 

Based on the results of the heat and mass transfer in a 
ventilated cavity, it is concluded that 

• Dl configuration: in cases where the source 
of contaminant is located near the heat 
source (CH3, CH6 and CH9), the lowest 
value of Taverage with respect to other 
positions of CO2 source can be obtained. 

• D2 configuration : for cases with a source of 
contaminant in the positions CHI, CH2 and 
CH3 (top) the best performance can be 
achieved from the standpoint of thermal 
and air quality. 

• From an air quality point of view, the 
Standard ASHRAE 62-2007 [13] states that 
the maximum concentration of carbon 
dioxide (CO2) inside a room environment is 
700 ppm, in order to be acceptable for air 
quality purposes. In general, the 
configurations Dl and D2 fulfill this 
requirement. 

In general, the air inflow is either near or far from the 
heat source and/or source of the pollutant, and the 
optimum location to situate the outlet flow of the 
mixture air-CCte is close to the heat source. 
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